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The evolution of two slow-rolling scalar fields with potentials of the form V = V(i4>~°' exp(— /3(/>™) 
is studied. Considering different values of the parameters a, jS and m, we derive several new 
inflationary solutions in which one fields just evolves in the background and is not important for 
the infiationary dynamics. In addition, we find new solutions where both fields are dynamically 
important during inflation. Moreover, we discuss the evolution of perturbations in both scalar fields 
and the spacetime metric, concentrating on the production of entropy perturbations. We find that 
for a large region in parameter space and initial conditions tensor modes are negligible and that 
adiabatic and isocurvature perturbations are essentially uncorrelated. 
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I. INTRODUCTION 

Since their advent just over twenty years ago, cosmological inflationary models have been developed to try to deal 
with several problems arising in the standard cosmology. The most important consequence of an inflationary stage in 
the very early universe is the development of perturbations in space-time and matter, which eventually evolve into 
the structures we observe in the universe today. The simplest model with an inflationary stage is where inflation is 
driven by a single scalar field, called the inflaton, with some potential V {(()). If the potential is flat enough, theory 
predicts that adiabatic perturbations are generated which obey the Gaussian statistics and have an almost scale-free 
spectrum (for recent reviews see e.g. and 0)- 

However, the most important drawback for inflationary cosmology is that there is no unique candidate for the 
inflaton field and its potential. Furthermore, according to our theories of particle physics, there are a large number 
of scalar fields which, potentially, could be important in the early universe. For example, hybrid infiation is a model, 
in which two scalar fields are important: whereas one field drives an inflationary epoch, the dynamics of the second 
field will end this period of inflation If infiation itself is driven by two or more scalar fields, the perturbations are 
no longer purely adiabatic or Gaussian 0]-C3- Moreover, the transition to the normal radiation-dominated epoch 
depends on the decay of the scalar fields into radiation and matter, which could infiuence the evolution of perturbations 
as well. Cosmological observations by the Wilkinson Microwave Anisotropy Probe (WMAP), the Planck Surveyor, 
2dF and SLOAN Digital Sky Survey will put strong constraints on the nature of the primordial perturbations and 
therefore put constraints on more complicated models of inflation and the process of reheating. 

Although formalisms have been presented in order to study two-field infiation (see e.g. ^^'^ ^3)j no 

explicit model has been used in order to examine the background evolution of both fields and the consequences 
for the perturbations. In this paper, we investigate the evolution of two scalar fields, each with a potential of the 
form V{(j)) = Vb0~" exp[— and examine the effect of the parameters a,(3,m on the evolution of the universe. 
This potential includes pure exponentials (a = 0), pure power law (/3 = 0) and combinations of the two. For 
a single scalar field, this potential was investigated in detail in and |l4j . Exponential potentials appear in 
Kaluza-Klein theories as well as in supergravity and superstring models (for a review, see e.g. ^^)■ Models based on 
dynamical supersymmetry breaking involving fermion condensates motivate inverse power potentials ; supergravity 
corrections to these models in turn predict potentials of the form above p^ . Furthermore, models motivated by 
higher-dimensional theories and/or scalar-tensor theories easily incorporate these potentials. For example, the form 
of the potential was used in [Tsj to examine inflation with a background dilaton field. 

The potentials just described were used both in inflationary cosmology as well as in models for dark energy, driving 
the observed acceleration of the universe at the present epoch (for a review, see Indeed, the potential above 

has interesting properties, allowing for scaling behavior and other attractor-like solutions (see e.g. and [IJ]). 
However, our goal in this paper is to obtain an understanding of the dynamics of two scalar flelds with this potential 
which drive an inflationary epoch in the early universe. This, in turn, allows us to study further the consequences for 
the perturbations, in the two scalar fields and the space-time metric, during the inflationary epoch. This can be seen 
as a first step to understanding the initial perturbations in the radiation dominated era, which eventually evolve into 
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the structures we observe today. We would like to emphasize, however, that in order to calculate the perturbations 
in the radiation dominated epoch, an understanding of the decay of both scalar fields is needed. There are some 
possibilities we would like to mention: 

• Both scalar fields decay, one into radiation and baryonic matter, the other field only into dark matter. The second 
field could, in principle, also decay only partly into dark matter, whereas the "rest" provides an explanation for 
dark energy. 

• Only one scalar field decays completely into matter and radiation, the other decays partly into some form of 
matter; the remains of the field play the role of dark energy. 

• One of the fields plays the role of the curvaton field [S^, i.e. it decays well after the inflationary epoch and is 
responsible for the curvature perturbation. 

There are more possibilities, of course, but it is clear that the subsequent evolution of perturbations in the metric and 
matter fields strongly depends on which of these possibilities is realized. The situation is now far more complicated 
when compared to inflationary models based on one scalar field alone. 

The paper is organized as follows: in Section 2 we present the equations of motion; section 3 discusses the case of two 
independent potentials for both fields, i.e. the total potential is given hyV = Vo^(j)~" exp[— /^^'"J+Vb^^X"" exp[— /Sx™]; 
in Section 4 we investigate a mixed potential of the form V = VoX~" exp[— In sections 3 and 4 we discuss all 
of the possible cases in detail. These two sections are necessarily technical. However, we have included a summary of 
our results at the end of each section. Readers who are not interested in the technical aspects may wish to refer to 
these summaries. The perturbations are discussed in Section 5; we present our conclusions in Section 6. 



II. THE FIELD DYNAMICS 



Our starting point shall be that of Einstein gravity with two scalar fields and a general potential V{(f>,x)- The 
action for this (in Planck units) is described by 



S = / d*xJ^c 



This action leads to the usual Friedmann equation and two Klein-Gordon equations for both scalar fields: 
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where we use the notation = dV/d(f> with its obvious extensions. 

We shall find it convenient to make the slow-roll approximation which is appropriate when the slow-roll parameters 
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satisfy |ei,77y| ^ — 4'iX- These five conditions are sufficient for us to generate an inflationary epoch where 

a{t) > 0. With this approximation the equations of motions (O reduce to 



= ^V{4>,x), 
3Hx = -V^. 



(4) 

(5) 
(6) 



For the most part, we shall see that slow-roll is an attractor in our model and that we shall be justified in making 
this approximation. 
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III. UNCOUPLED POTENTIALS 

In this section we begin our investigation by considering two uncoupled potentials of the form 

V{cj,,x) = Vo,r"' exp + Vo^x'"' exp [-/32X"1 ■ 

In this case the slow-roll parameters read 



(7) 
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(8) 



and similarly for with rj^-^ = 0. It is clear that, with a judicious choice of initial conditions, these can always be 
made small. It is also clear that = 1 will be a transitional value because, for > 1, the slow-roll parameters will 
grow. We now begin to examine the evolution of the fields in this regime. We shall assume that (3i > but make no 
further assumptions about the other parameters in the potential terms. 



Large Field Values 



A. mi, 1712 > 

With TOi,m2 > 0, the exponential part of the potential, for large enough (j),x at least, will always dominate over 
the power law. Dividing equations H5I6|) we generate the equation 



(9) 



We are unable to integrate this directly but are able to do so to leading order in the fields by means of integration by 
parts. This then generates for large (f), x 



2^2 ,2mi-2T 



yo>~"^ exp [-/3i</.™^] = f3lmlx^"^^-^Vo^x-''' exp [-/32x" 



(10) 



One then takes natural logarithms to find a simple expression between the fields. Since > the terms proportional 
to the powers of the fields will eventually dominate. This gives us 



to leading order. Then using and we are able to show 



(11) 
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Substituting back into (jSJ we must solve the differential equation for 
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(12) 



(13) 



(14) 



to leading order in (/). For large </>, we are able to integrate this expression under two different regimes, depending 
on whether mi < m2 or mi = m2. This corresponds to one field giving the dominant contribution to the Hubble 
parameter and both giving equal contributions. This is clear from the denominator in (|14|l . We shall take both of 
these cases in turn. 
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1. mi < 1712 

In this case it is easy to see that the (p field is the dominant contributor. Solving for (j) we integrate by parts 
to show 



_3_ 



-2mi+2 



exp 
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Immediately we are able to see that large (j) implies large t and so inflation will occur at late times with this set up. 
We are unable to invert this directly. However by taking In, we see that the exponential quickly dominates, and so to 
leading order 
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Writing H15() in the form 



= In 



+2mi-2 



(17) 



we are able to substitute back into \i<6\ to obtain a solution for (p{t) to next order in t. A little thought reveals why 
we want and need to do this. When we come to evaluate the Hubble parameter as a function of time we need to take 
the exponent of cf). Every time we do this we must take the next to leading order expression. This is given by 
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Then we plug this into (|12|l to solve for H{t) and then integrate to get the scalefactor which gives us, 
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(19) 
(20) 



In the evolution of the fields we may drop most of the coefficients as they simply result in constant terms in the 
solution. This gives the final solution 



(21) 
(22) 



Ht) = ( ^Int 



Xit) = [j^lnt 



Now because < 1, we see that the scale factor grows faster than power law. This regime is effectively equivalent to 
a single scalar field driving the dynamics of the universe as the second scalar x sits in the background. Unsurprisingly, 
the solutions are in direct agreement with section III. A of [T^. 



mi > 7712 



This is equivalent to the previous section by symmetry. 
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mi — 7112 = rn 



We now consider the case where both the fields give equal contributions to the evolution of our universe. Since, the 
second field is now significant, we expect to see something new. This case includes that of two straight exponential 
potentials which has already been studied in We are able to integrate ()14|l to leading order to show that 
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Once more the evolution of ^(0) demonstrates that inflation occurs at late times for large 0, x- We then proceed in 
the same way and the solutions we get are as follows: 
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It should be noted that in the two exponential potential case this reduces to power law inflation 



ait) - t 



(24) 
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which is agreement with the dynamical system analysis of 23]. The "new" feature is that we generate more inflation 



than we would get with one scalar because of the extra factor, exp 
This concludes the analysis for the cases mi, TO2 > 0. 
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in the solution for the scale factor. 



B. mi < 0,m2 > 

It is possible when to^ < and is also an even number, to generate a minimum of the potential at the origin. The 
nature of the slow-roll parameters dictates that this will certainly give a finite amount of inflation and that we may 
not get slow-roll here. We shall assume, initially at least, that the field begins in the part dominated by the power 
law part of the potential. 



1. Qi/0, -4 

Following the same procedure as the previous section we are able to deduce that 
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ai(p 2 exp 



Taking the 4> exponential as approximately constant, we are able to solve this to give. 
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Solving for the remaining fields and parameters we deduce 
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In this instance 



fi{t) = exp 



1+4 




(35) 



and rapidly tends to unity as time increases. With this set of parameters, studying H30() one realizes that it is possible 
for inflation to occur in a different manner. It is helpful to write the scalefactor in the form 



a{t) cx exp [ct^ fi{t)\ 



(36) 



ai > 0: In this case, ai(ai + 4) > and so it follows, from equation 1301 that t((/)) ^ oo as increases. Therefore, we 

must generate inflation at late times. In addition, c>0,0<fc<l giving us a relatively slow expansion rate. 

— 4 < ai < 0: In this case ai(ai + 4) < and we get c > 0, 1 < /c < oo. As ^ oo, t[(j)) —oo and so inflation now 

occurs at early times. Potentially this can also give us a lot of inflation as the scale factor grows faster than e*. 

ai < 0: Once more we get ai{ai + 4) > 0, however we now have a negative power of cj) in (|3()|l and so t{(j)) — > as 

increases. This generates inflation at early times once more. In addition we have the constraints c, /c > 0. 

Once more the setup is equivalent to that of a single scalar as the second scalar makes minimal contribution to the 
universal dynamics. We generate identical results to section III.B of [T3| . 

Let us now return to the case where x dominates. This gives a Hubble parameter of 
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Integrating up, again regarding the exponential as almost constant we get 
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We are not able to invert this equation in order to obtain a meaningful description of (^{t). 



2. Qi = -4 

The evolution equation for (j) H29(l still holds - provided, of course, that we are still in a dominated regime. 

(i = -4 
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Then integrating we retrieve 
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Notice that we generate a In rather than a power of (j) in the solution. This propagates through the calculation and 
will generate qualitatively different behavior in the Hubble parameter and scalefactor. However, we still see that 
t{(f)) — *■ — oo for increasing </) and so it is an early effect. 



exp 



m 

xit) 

H{t) 
a(t) oc exp 



P2 V 3 
Vn 



exp 



1 

-— exp 



where 



hit) = exp 



3 

}\ 
3 



Pi 

— exp 



t exp 

2 



-4toi 



Vn 



hit) 



-4mi(^) t 



(42) 

(43) 
(44) 
(45) 

(46) 



It is clear that fsit) 1 rapidly as t increases. Looking at equation H41|) . it is clear that t decreases for increasing 
so the inflation we generate is an early time phenomenon. 



3. Qi = 



In this special case the (j) field will quickly dominate once more, as the potential for x decays quickly. Thus we get 
behavior very similar to that for one field, with x just in the background as in |l4l |. Solving in the usual manner gives 
us 
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Inflation here occurs at early times. 
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C. mi, m2 < 
Proceeding in the usual manner we are able to show that 
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(53) 



It is then clear, for example, that —2 < ai < 0,2 corresponds to cf) domination. We shall consider this case and that 
of equal contribution since x domination is identical in behavior to the former. 



1. Q2 > Qi > — 2; a2 < — 2,ai > 02 



If is clear that the (f) field is dominant here and so the results from section (jlll B 1\ hold for 0(i), H(t) and a{t). In 
this case we do get different evolution for x(t) though, given by 
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2. ai = a2 = a 

This is the only case where we do not get domination by one field over the other. Again the results will be very 
similar to that of section (|111 B l|l except that we will generate more infiation because the Hubble parameter will be 
larger than before. Then we find that: 
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where 



/4(i) = exp 



(60) 



fiit) —>■ 1 rapidly as t increases. This gives us inflation at late times once more. 



Small Field Values 



A moments thought reveals that in order to generate slow-roll at small values of either x or cj) we must have 
ai — 0,m> 1. This can be seen easily by examining the slow-roll parameters Q. Our analysis will then follow a 
similar pattern to before. Again, we shall see that one field will become dominant if jtii ^ m2 and we proceed as 
before. 



D. mi < m2,mi ^ 2 
Once more the will provide the dominant contribution and so 



exp 



Then we find that 
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This is identical to the case we considered previously in section pil B 3|1 and so all the results there still hold. The 
one difference is that the evolution of x changes slightly. 
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One should note that as ^ 0, t decreases and so we are looking at an early time effect. This is equivalent to section 
IV. A of 01 . The cases where either one or both of rrii = 2 are easily covered and does not result in any analytical 
difficulties. The behavior is not markedly different from the above. 



E. 2 = mi < 1712 
The (j) is still dominant and so we can integrate (|52|l to achieve 
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Solving the equation for small we find that 
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Note that although the behavior x{t) ~ {—Kit)'^^,Ki > appears problematic, this would be resolved by including 
the integration constants. This would then give us 
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3 ; 1^0^/32^2(2^7712) 

Now since 4'{t) in a monotonic increasing function, the slow-roll conditions will be violated within finite time and 

'2- 



so our solution will no longer be valid. This should occur before ^ = -j ("vf^) Y!hL f^^'^'^i'^ 



since there are no 



physical grounds upon which to rule out this solution. 



In this instance we find that 



It is not too difhcult to show that 



F. mi = 7712 = 2 
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V3 



(68) 



(69) 



It is clear that as — > we have t — oo and so inflation will occur at early times. Then, proceeding in the familiar 
manner, we can show 
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where we have treated all the exponentials as being approximately equal to unity. 



(70) 

(71) 

(72) 
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Leirge and small field values 



If we take one of the fields large and one small we will always have the large field dominant and so this relates to 
the cases considered already. For example, small x means that 7712 > 1 and so 7711 < 7772. Therefore the (j) held will be 
dominant. This applies to all possible cases. 



G. Summary for Uncoupled Potential 



We have seen that for two fields with uncoupled potentials often one field will be dominant, and we return to an 
effective single field theory - for the background at least- as studied in There are, however, several cases where 
the second field is important and we have demonstrated them. We have seen that it is possible to generate infiation at 
both late and early times and for large and small values of both of the fields. As a result, these background solutions 
do not differ qualitatively from the single field set-up |l4j . The results are summarized in the table below. If mi < 777,2 
then the (j) field is dominant so we only outline the results for mi. If mi = = m the same qualitative behavior 
exists except that we generate more inflation than we would normally with one field. 



11 



TABLE I: Summary of Inflationary Behavior for uncoupled potentials 



Early Times 


Late Times 


All Times 


mi < 0, ai < 


< mi < 1 


mi = 1, ai = 


or 


or 


or 


mi > 1, ai = 


mi < 0, Qi > 


mi = Qi = 



IV. A COUPLED POTENTIAL 



We now turn our attention to a potential of the form 



(74) 



and repeat the processes of the previous section. The motivation for such a potential arises, for example, in higher- 
dimensional theories of the form 



d^x\/—gN 



(75) 



Writing the A^~dimensional metric in the form 

dsjf ^ gl^^dx^dx" + b^dsl, (76) 

and compactifying the theory on a torus which has zero curvature, one ends up with a potential of the form H74|l . 
where x esc In 6 describing the size of the extra dimensions (see e.g. 0). 

We still have the three coupled equations 14lt)|) to solve which are equivalent to 
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There are certain cases where it is not clear that equation (I80|l means anything- for example m > 2, a > 0- since 
this would generate an imaginary value for at least one of the fields. At this stage one should remember that we have 
omitted all integration constants up to this point. As there are no physical grounds on which to rule out such cases, 
when one puts back these constants (dependent on the initial conditions), any such problems may be overcome. 

Before we begin, we shall find it instructive to examine some plots of the potential in question. These are shown in 
figure 1 where we have plotted V[(f),x) = 4'^^ exp [— 2x^°'^] respectively. 

Large field Values 



A. 1 > m > 

This guarantees us slow-roll for large values of the fields. We generate the following equation for x 

-/3X" 



X = ( y I Pm[ ^- I X * X exp 

We are then able to solve this to produce the following expression, 
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FIG. 1: Plot of the Coupled potential. 
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Proceeding in a similar, way we are able to solve the system producing, to leading order in t, the results below. 
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(83) 
(84) 
(85) 
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One should note that this gives very similar behavior to section Ijlll A ip . This should not be entirely unexpected 
because the exponential part of the potential will always dominate for m > 0. 



B. m < 

Due to the nature of the slow-roll parameters, we must once more consider large values for both of the fields. Now 
solving for 0, treating the exponential approximately constant, we find 



Substituting for either field using l|Sn|) we able to invert this finding 
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One should also note that the onset of inflation occurs in analogy with section B 111 . It is also worth pointing out 
that this gives the solution one would expect when to ^ 0. That is, we generate the same results as (IIII B 1|) except 
that the x held freezes out. This should not surprise us since the potential becomes equivalent to that of an inverse 
power. We then flnd the form for the Hubble parameter 
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13 //3m(2- m)y- 

y \ 2^ 



-(a + 4)i 



(2-r,i)(Q+4) 



(90) 



where, of course, the exponential part tends rapidly to 1. Solving for the scale factor we deduce to leading order 

4 + a fVn\^ ( fVo\^ a 
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where 



giit) = exp 
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It is obvious that the special case of a = gives us just a single scalar field that evolves. This has already been 
studied. We must however consider what happens when a ~ —4. We generate the following equation 



1 / 3 



4 \Vn 



[In 0] exp 



Then solving for the fields we find 

m 



where 



exp 



exp 



X(i) 

m 

a{t) oc exp 



2 - m V 3 



exp 



, 13 

t exp 

2 



8m 



2-771 V 3 



t 



1 

-- exp 



92{t) 



g2{t) = exp 



2 exp 



2 - 777 V 3 



(93) 

(94) 
(95) 
(96) 
(97) 

(98) 



which rapidly tends to unity. 



C. 1 < m < 2 



If one considers the slow-roll parameters we may expect only to generate inflation for small values of x- However, 
we must still take large <j} for slow-roll to be a valid approximation. In this instance, we are able to approximate the 
exponential by unity. Provided m ^ 2, we generate the same set of results as for 777 < with large x, as in section 

In the special case where 777 = 2, we get slightly different behavior. We summarize the results below: 



where 



x{t) 



a 



{a + 4)t 



exp 



a(t) oc exp 



exp 



--exp 



m(3 
a 



-(a + 4)t 



(99) 
(100) 

(101) 
(102) 



gait) = exp 



f3 



■ exp 



777/3 



Vo 



(103) 
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D. Summary for the Coupled Potential 

The summary of the inflationary behavior is efltectively equivalent to Table ^ with the obvious changes to the 
parameters. We are able to generate inflation at late, early or all times by tuning the parameters in the potential. If 
m > 0, the exponential part is dominant and the solutions produced are similar to that of a single scalar, x with the 
second field, (/>, in the background. When to < 0, the exponential part is subdominant and the results resemble those 
for a single scalar, 0, with a potential. The quantitative modifications are small and quickly decrease with time. 

In all cases one of the fields provides the dominant contribution to H{t) and a{t). The effects of the second field 
are present but decay quickly so that the field is only in the background. 



V. COSMOLOGICAL PERTURBATIONS 



We now turn our attention to the evolution for cosmological perturbations in both of the scalar fields and the 
spacetime metric during a period of inflation. Again, we assume that the fields are slow-rolling, implying that the 
parameters © are small. Furthermore, we shall assume that both fields contribute approximately equally to the 
expansion rate. Only then can we expect new features in the perturbation spectra resulting from the second field. 

A convenient formalism to study both adiabatic and isocurvature perturbations in inflation was presented in 
which we will use here. Instead of working with the fields (j) and x it is useful to perform a rotation as follows: 

5(7 = (cose')50+ (sin6l)(5x (104) 
5s = -(sin6i)(5<?!)+ (cos6i)(5x, (105) 



sin 9 = (106) 



with 



a is called the adiabatic field and s is called the entropy field. The motivation for their names becomes clear when 
one considers fluctuations of them. 

The line-element for arbitrary scalar perturbations of the Robertson- Walker metric for a spatially flat universe 
reads, (using the notation of 8]) 

+ a^[{l-2^P)5.,j +2E^,j]dx'dx^ . (107) 
The gauge-invariant curvature perturbation, defined as 

7^ = ^ + S, (108) 

p 

is, on very large scales, constant for purely adiabatic perturbations (see e.g. |25j| and references therein). However, 
entropy perturbations are a source for the curvature perturbation (|108|l . The entropy perturbation between two 
species A and B, defined as 

5=^-^, (109) 
riA riB 

where are the number densities of the particle species i can evolve in time, even on superhorizon scales. Therefore, 
it was argued that on very large scales in general we have the following equations describing the evolution of TZ 
and S: 

n ^ -fHS, (110) 
S = 5HS. (Ill) 



For the case of the two slow-rolling scalar fields and in the spatial flat gauge {ip — 0), TZ and S are given by 
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and 



H 



S=^. L = ^. (113) 

For S we use the same normalization as in 12]. Fluctuations in the field a are adiabatic perturbations, whereas 
fluctuations in s are entropicperturbations. On very large scales (fc <C aH) and in flat gauge, the evolution of 
fluctuations are described as Q 

(Sa)- + 3H{Say + (v^^ - 6^^ Sa = -2V„A + + 2(0(5s)- - 2^%s (114) 

and 

(5s)- + ZE(8sy + [Vss - 9^^ Ss ^ -2^ [cr((5cr)' " " ^^^a] . (115) 
The metric perturbation A can be obtained from Einstein's equation and is given, in the flat gauge, by 

HA^ATrG(4>d(l> + xSx) ■ (116) 



An important point is that 9 must be non-vanishing in order for the adiabatic field to be sourced by the entropy field. 
If 9 is constant, the entropy field does not contribute to perturbations in the gravitational potential. 

In 12'] the formalism presented was applied to slow-roll infiation with two scalar fields and it was shown that 7 
and 5 in (|110|l and (|lll|l are given in terms of the slow-roll parameters 

7 = -2?fos, (117) 
6 ^ -2e + 7]„„-r]ss, (118) 

and are, therefore, specified by the potential V{4',x)- Iii the last two equations, the slow-roll parameters are con- 
structed from the usual slow-roll parameters for (f> and x © s-iid ^-re given by 

'^liy) + (119) 

and 

V<j<y — V<P4> cos^ 9 + 2ri^x cos 9 sin 9 + rj^^ sin^ 9, 

Vss — rj^cf, sm^ 9 — 2r]^^ cos 9 sin 9 + rj)^^ cos^ 9, (120) 
ViTs = (Vxx ~ V4>4>) sin 6* cos 6* + rj^^ {cos^ 9 ~ sin^ 9) . 

The time evolution of TZ and S between horizon crossing and some later time is given by 



n\ 1 Tns] n 



(121) 

where the asterisk marks the time of horizon crossing. The transfer functions T-jis and Tg^ are given by 

Tss{t,t.) = exp(^^ 5{t')H{t')dt'^ 

Tns{t,U) - / i{t')Tss{U,t')H{t')dt'. (122) 



Thus, in the case of two slow-rolling scalar fields, the transfer functions are completely specified by the potential 
through the slow-roll parameters. 

We are now in a position to calculate some of these transfer functions for our general potentials. We shall consider 
the most simple cases only in order to show how the background solutions derived in Section 3 and Section 4 can be 
used in order to study perturbations. 
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A. Analytic Calculations 

We begin with analytic calculations for two simple cases. The first case considers two exponential potentials (a 
special case of assisted inflation). The second case is a simple model of a coupled potential. 



1. Two Exponentials 

Let us now take 

X) = V^^e-P^* + V^^e-P-^. (123) 

Although this setup was studied earlier pTj . we would like to check the formalism for the perturbations. We return 
to the first approximation Hl()|l to make the calculation. Then in this instance we find 



0- 

P2 



Vo,^-^''^ = ^^o.e-/'^^ (124) 



^ X- (125) 
Pi 



These two equations enable us to calculate all the parameters. We find that 

02 . „ 



cos9= , sin6> = . (126) 



This shows that we get a constant angle in the phase plane. A quick glance at equations 1)1141115(1 reveals that the 
adiabatic and entropy perturbations decouple in this case. We should expect our analysis to reflect this. It is also 
straightforward to show that 

2e = ... = .. = ^, (127) 

Vas = 0. (128) 

The functions "f{t), S{t) then turn out to be constant given by 

7 = 0, (129) 
* = -^<0- (130) 

Then we simply compute the transfer functions using H122|) to see that 

Tns{t,U) = 0. (132) 

Now since 5 < 0, it is immediately obvious that Tss{t, t*) ^ as time passes, provided t > t*. Secondly, T-jis{t, t^) = 
0. This means that TZ remains constant and 5 — s- from H121|) . Therefore, in this case we only expect adiabatic 
perturbations and no entropy perturbations if inflation lasts long enough. This is in agreement with the results in 
[27| . One should add that the transfer functions are calculated from the time of horizon crossing. However, when we 
calculate the transfer parameters, ^{t), 6{t), we assume that k/aH <^ 1. This is only true a few e-folds after horizon 
crossing. This means that our calculated transfer functions will be accurate a few e-folds after horizon crossing and 
so, in effect, they have the wrong initial conditions. Whilst this does not affect their predictions, this means they 
differ from what we observe numerically. In this instance there is some evolution of TZ before it freezes in and some 
evolution of S before if follows the decay predicted. This can be seen in Figure [21 in the first few e-folds after horizon 
exit. 
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2. A coupled example 



The next simplest case will prove to be that with a coupled potential with a straight exponential. That is 

^^or^e-^'^. (133) 
We shall find it most convenient to work in terms of the field (p before substituting in its solution (|83|) to compute 



our integrals. We find that our 9 now evolves according to 

a 



cos 9 

We also have the slow-roll parameters 

V<t>4' - 

It is then straightforward to show that 



:, sin 9 



v/a2 + /3202 



a[a + 1) 



a2 /?2 



75 \-2(j + —0 



(3 



a2 + /3202 Q, 

We then find that the transfer parameters are given as 

2/3^2 



q;2 + /3202 0- 



lit) 

S{t) 



(a2 +/3202) 



1 



q;2 + /3202 



02 Q/ 



(In i) 2 (a + 4 hit) 
_ /32(Q-41nt) 
~ 4(lnt)(a + 41nt)' 



(134) 

(135) 
(136) 

(137) 
(138) 

(139) 
(140) 



We have also verified these relations numerically. Now it remains to compute the transfer functions which we are able 
to do exactly. From (|85|) we see that 



This then gives 



where 



,* (a-41nt') 1 
~ 2(lni')(« + 41nf)t^ 



Ini 



Int* 



In 



a + 41nt 



a + 4 In 



Plugging this back in we find that 



Tssit, U) 



Then using this result 



Tns{t,U) = (3^ a--'- 



' Int ■ 


2 


'a + 41nt*' 




Ini* 




a + 41nt _ 




+ 41nt. 


/■* 2 


1 


(lni») 


1 
2 


it. /32i' {a 


+ 41nt')^ 



2a 5 



a + 41nt, 



dt' 



a2 / 1 



1 1 



4 a + 4 In t' 



a + 41nt, 1 



X\nU)^ (lni,)2 a + 



4 hit 



(141) 

(142) 

(143) 
(144) 

(145) 

(146) 
(147) 
(148) 
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Thus at late times we see 

a 2 1 

Tss ^ 0, Tns -> -X r- (149) 

2 (Int.) 5 

Again, this imphes that at late times TZ remains constant, although its value on a certain length scale depends on the 
time of horizon crossing . 



B. Numerical Calculations 



In principle, with the background solutions given in Sections 3 and 4, it is possible to calculate any transfer function 
that we wish to consider. However, in practice, this requires making a large number of assumptions in order to do 
the integrals analytically. Therefore, to get precise answers, we will now resort to a numerical investigation of the 
perturbations. We begin with the case of two exponential potentials, where we can check our numerical results with 
analytical solutions. Then we study two inverse power potentials, then the case of a mixed potential and finally the 
coupled potential. The initial conditions for background are chosen such that both fields contribute approximately 
equally to the expansion rate. With this prescription a change in the initial field values is equivalent to changing the 
amount of inflation with the same parameters. 

We choose to run our simulations in terms of the variables Sa and Ss - this gives a more accurate solution. The 
equations for their evolution in the spatially flat gauge are given by 



Sa + 3HScT 



3HSs 




(150) 



(151) 



The general recipe is as follows: we evolve the background through 100 e-folds of inflation to make sure that the 
fields are in the slow-roll regime. We then switch on the perturbations so that each mode starts inside the horizon 
with k — lOOOaif at different times. We then follow the perturbations for another 60 e-folds after horizon exit, where 
inflation is assumed to end. 

To set the initial conditions one treats the modes, Sd and 6x, as independent stochastic variables deep inside the 
horizon. With this prescription, it can be shown that 0,0] 



aScf) ■ 



1 



-ikr 



2k 



(152) 



deep inside the horizon, and similarly for Sx. In the last expression, dr = dt/a is the conformal time. To calculate the 
spectra numerically, we use the method described in |l(Tl |. One makes two runs: the first run begins in the Bunch-Davis 
vacuum for dcj), Sx = Sx = 0, and the second in the Bunch-Davis vacuum for 5x, S(j) = Scj) = 0. The power spectra 
are then given by 



2^ 



1,3 



Vc 



k^ 

2^ 



|7^l5l + 7^252| 



(153) 
(154) 
(155) 



where the subscripts 1, 2 are the results for the two different runs. Furthermore, for the tensor modes one finds 

2 / X 2 

nr - - I - I ■ (156) 



We also define 
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Note that our normalization of rr is 16 times smaller than that used in '2^. As a test for our code, we are able to 
check it against the results for assisted inflation with two exponential potentials W\ . The analytic results give |23 



-, ^A^lH2 ^^KQ\ 



^{^1 + Pi - 3/32/31) (/32 + f32 _ 1^2^; 



Finally, in the case of two-field slow-roll inflation there are two consistency relations: 

rr = -^{l-rD, (160) 

{nc-ns)rr = — — {2nc - nu - ns) ■ (161) 

In the models with runaway potentials considered here we generally find that the spectral index of the correlation 
spectrum is much larger than the indices for the adiabatic and isocurvature perturbations. Furthermore, the pertur- 
bations are essentially uncorrelated so that the two consistency relations are effectively equivalent. We find that they 
are both satisfied to within a few percent either side. 

1. Two Exponential Potentials 

Let us discuss first the case of two exponential potentials. We have seen from the analytical treatment that in this 
case the entropic perturbations decay. We now give the fields a power in the exponent, Vicj), x) = Vq^ exp [— /3i0™i] + 
Vo^ exp [— /32X™^]. Using equation l(Tn|. one can show that 

tan..i = ^H!!^. (162) 

Then, provided the two fields have different potentials and initial conditions, one would expect the entropy pertur- 
bations to be sourced because 9^0. To gain further understanding, if one makes the field redefinition, ip — x™^ we 
may write the Lagrangian as 

= i^^^ (9^)2 + exp [-/32^] . (163) 

7712 

We see that we generate a non-canonical kinetic term. It has been shown, , that such a term sources perturbations 
in 6s and hence in entropy. Because these source terms are proportional to derivatives of this non-canonical factor, 
these decay for 771.^ < 1. Furthermore, since the potential is then reduced to that of straight exponential, we would 
still expect the entropy perturbations to decay, the difference being its duration. 

Let us now consider some examples and begin with two exponentials. We set the normalization such that TZ ^ 10^^ 
at the end of inflation. 

In our first example, we set /3i = 0.05, /32 = 0.10, 7771 = 7^2 = l,Vo = 10^^°. Figure [21 shows the evolution of 
TZ and S for a given mode from the time of horizon crossing. We can see that the entropy perturbation begins to 
decay almost immediately as one would expect from our previous calculation of the transfer functions. Further, one 
can see that there is significant evolution in the first few e-folds after horizon crossing, N^, « 8. Following this TZ 
remains constant as expected. As a second example we consider a case with a power in one of the exponentials. The 
parameters are as before setting 77ii = 0.3. Above we have shown that this should source the entropy perturbation 
which is realized in Figure |21 Let us now observe how the spectral indices behave for all of parameter space. We set 
/3i = 0.05, /32 = 0.10 and vary the powers 7711,7772. We maintain the initial conditions (/>o = l,Xo = 1- The results are 
summarized in Figure |3| We find that nc always runs significantly so we do not include it here. 

We observe in Figure O that the spectrum for the tensor modes is flat- to the extent that is indistinguishable from 
a flat spectrum with present observations- and the amount of tensor production is very small. Again this is below the 
level of current experiments. In general, the spectral indices are functions of the slow-roll parameters which decrease 
in time with run-away potentials. Thus the more inflation we allow the smaller the indices will be and the closer 
to scale invariance we will be. In fact, in all cases we can get arbitrarily close to a Harrison-Zel'dovich spectrum 
by allowing an arbitrarily large amount of inflation. Furthermore, one may expect the plots to be symmetric but 
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FIG. 2: Evolution of TZ and S for a wavelength that leaves the horizon, A^, = 8, 60 e-folds before the end of inflation. The 
left figure shows the case /3i = 0.05, P2 = 0.10, mi = m2 — l,Vo = 10~^° with initial conditions (t>o = l,xo = 1- On the right 
hand side we now set mi = 0.3 leaving the other parameters unchanged. The plots show |-Pr| 2 and jPsI ^ . 



we remind the reader that /3i 7^ (32 so the two fields always have different potentials. Our results do not completely 
agree with the analytic results given above. The reason is that we have not yet reached the attractor solution for 
two exponential potentials. However, we have verified that once we are on the attractor, the analytic and numerical 
results agree sufiiciently. 



2. Two Inverse Power Potentials 



To gain further insight here, it is helpful to use the field re-definition, e = x Then, 

C^ = ^ei^{d4>f + V^^e-''. (164) 

Therefore, this case corresponds to two scalar fields with non-canonical kinetic term and exponential potential. Be- 
cause the non-canonical kinetic terms are a source of entropy perturbations, we expect that the entropy perturbations 
to be more strongly sourced when compared with the previous example of two straight exponentials. The data are 
shown in Figure^ Compared with the case of two exponential potentials, the contribution of the entropy to the 
total perturbation is larger, but we see that the entropy modes are always smaller in magnitude than the curvature 
perturbation. Furthermore, adiabatic and isocurvature perturbations are effectively uncorrelated. The cleft in the 
plot arises from the background fields having the same evolution as they each have the same a. It is the clear that 
from our definition of Pc that this quantity should be almost zero in this cleft. Once more the spectral index of the 
tenors is very small. 



3. Mixed Case 



We now consider the case of one field having an inverse power potential and the other an exponential potential. 
Again, we can transform the inverse power law potential into a exponential potential. It is then clear from the 
discussion of the last two subsections that the results for the spectra must be similar, as one can observe in Figure|Sl 



4. Coupled Potential 



Once more the features observed here are very similar to the cases studied so far. We set (3 — 0.05 and vary the 
other parameters. The data are shown in Figures IHI 
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-0.05 
-0.1 
-0.15 





FIG. 3: Properties of the spectra for two exponentials. We set /3i = 0.05, /32 = 0.10 and vary the powers m\,m2 and we 
maintain the initial conditions (f>o = l,xo = 1- AH quantities are defined in the text. We remind the reader that the scale of 
Ps is in general arbitrary but here it is equal to Piz at horizon crossing. 



23 






FIG. 4; Properties of the spectra for two inverse power potentials. We vary the powers ai,a2 and we maintain the initial 
conditions (f)Q = 1,xq = 1- AU quantities are defined in the text. We remind the reader that the scale of Ps is in general 
arbitrary but here it is equal to Pr at horizon crossing. 
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FIG. 5: Properties of the spectra for the mixed potentiaL We set /3i = 0.05 and vary the powers mi,a2 . We maintain the 
initial conditions (jio = 1, Xo = 1- AH quantities are defined in the text. We remind the reader that the scale of Ps is in general 
arbitrary but here it is equal to Pr at horizon crossing. 





FIG. 6: Properties of the spectra for the coupled potentiaL We set /3 — 0.05, vary the powers m, a and we maintain the initial 
conditions 4>q = 1,xq = 1. AH quantities are defined in the text. We remind the reader that the scale of Ps is in general 
arbitrary but here it is equal to Pr at horizon crossing. 
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5. Compatibility with WMAP 

In general, all our models generate almost scale-invariant spectra with minimal tensor production. This also allows 
both red and blue spectra as seen in Figure Furthermore, we are able to get arbitrarily close to scale-invariance 
by allowing more inflation. In addition, we find that the perturbations produced are essentially uncorrelated. 

Let us comment on the compatibility with current experiments. The WMAP results have provided strong constraints 
on the initial power spectrum and hence on inflationary models |2^. If one examines Figure[7|one sees that it possible 
to have 

0.9 <nn< 1.3 at 2a. (165) 
We include the equivalent plot for our models. 




1.014 



n_R 

FIG. 7: Allowed region for curvature index from the WMAP data alone. We scale our tensor/scalar ratio to agree with the 
WMAP team's plot. One can clearly see the degeneracy in the predictions of rr and n-n among the possible models. For example 
as m — > in the coupled case we approach the results of the inverse power potential with the same dominant a. This is what 
we would have expected from the background solutions of sections 3 and 4. 

One open issue we have not discussed is the end inflation. This may be caused by another scalar field. This process 
is, in general, non-adiabatic and would generate further entropy production which would add to that from the two 
fields. Hence, the spectra deep in the radiation dominated epoch depend on the details of reheating. Without a 
detailed study of this mechanism it is impossible to draw firm conclusions from the entropic values in all cases. 

VI. CONCLUSIONS 

In this paper we have studied in considerable depth inflationary solutions generated by two scalar fields and 
potentials of the form V = Vb0~" exp(— /30™). This potential is very general and encompasses potentials motivated 
by supergravity and string theory and more general theories motivated by higher-dimensional gravity. 

We have seen that it is possible to generate a large range of different inflating universes by varying the parameters 
in the theory. We find that slow-roll is a valid approximation by means of numerical calculations and we are able to 
verify that the analytic solutions presented in this paper hold. For two uncoupled fields (Section 3), we found that, 
often, only one of the fields is dominant. However, we also found solutions where both fields are important for the 
infiationary dynamics. For two coupled fields with the above potential (Section 4), we have found that only one of 
the fields dominates inflation, whereas the other field is only a background field. This case might be interesting if the 
background field acted as a curvaton field, for example. 

We have used some of the background solutions we found to describe the evolution of the perturbations in two 
specific cases and, potentially, one would be able to find the transfer functions for all the solutions discussed here. 
We have demonstrated this with two examples. In general, however, one has to resort to numerical methods in order 



27 



to evaluate the transfer functions. The solutions presented here thus allow one to make predictions about the relative 
size of TZ and S for a large class of inflationary scenarios. 

We made use of numerical calculations to investigate consequences for cosmological perturbations. We discussed the 
case where both fields contribute approximately equally to the expansion rate during a period of slow-roll inflation. We 
found that our models produce scale-invariant spectra with a small scalar-tensor ratio. By varying the parameters of 
the potential, it is possible to produce both red and blue spectra for the curvature perturbation. This is still consistent 
with observational constraints. Although the exact details depend on the initial condition of the fields, we don not 
expect our conclusions will be changed qualitatively. In the case of double inflation, one can generate uncorrelated 
[^,|^ as well as correlated perturbations, |29j |. In contrast, in multi-field inflationary models with run-away potentials, 
the perturbations are always uncorrelated. 

In order to make predictions for the anisotropics in the cosmic microwave background radiation, one has to follow 
the perturbations into the radiation dominated era, which involves a detailed calculation of the decay of both fields. 
Ending inflation in the class of models discussed here will involve presumably a third field, which becomes important 
just at the end of the inflationary stage. 

In future work we will use some of these insights gained in this work in order study inflationary models in models 
with extra dimensions |30|. 
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